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Escape fraction of ionizing photons from galaxies at z = 0—6 



O 
O 



o 



> 
(N 

in 
in 
o 

O 

Oh 
6 

o3 



13 



Akio K. Inoue 1 *, Ikuru Iwata 2 , and Jean-Michel Deharveng 3 

1 College of General Education, Osaka Sangyo University, 3-1-1, Nakagaito, Daito, Osaka 574-8530, Japan 

2 Okayama Astrophysical Observatory, National Astronomical Observatory of Japan, Kamogata, Okayama 719-0232, Japan 
3 Laboratoire d'Astrophysique de Marseille, Traverse du Siphon, BP 8, 13376 Marseille, CEDEX 12, France 



Accepted Received; in original form 



ABSTRACT 

The escape fraction of ionizing photons from galaxies is a crucial quantity controlling 
the cosmic ionizing background radiation and the reionization. Various estimates of 
this parameter can be obtained in the redshift range, z = 0-6, either from direct 
observations or from the observed ionizing background intensities. We compare them 
homogeneously in terms of the observed flux density ratio of ionizing (~ 900 A rest- 
frame) to non-ionizing ultraviolet (~ 1500 A rest-frame) corrected for the intergalactic 
absorption. The escape fraction is found to increase by an order of magnitude, from a 
value less than 0.01 at z < 1 to about 0.1 at z > 4. 
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1 INTRODUCTION 

Although recent observ ations have revea l ed the outline o f 
the cosmic reionization (Pag e et alJl2006r iFan et alj |2006). 
its detailed history and the nature of ionizing sources are not 
yet fully understood. Early forming galaxies can be strong 
ionizing sources, depending on the escape fraction of their 
ionizing photons (/ CS c). This / esc is therefore a key quantity 
for understanding the cosmic reionization process but its 
typical (or effectiv e) value is not yet clearly established (e.g., 
Ilnoue et al]l2005T) . 

The goal of this Letter is to put together all existing 
information on the amount of ionizing photons released by 
galaxies into the intergalactic medium (IGM), both from 
direct observations and from indirect estimations based on 
the effects of the IGM ionization itself. In order to make an 
homogeneous comparison, we introduce a quantity derived 
from / csc , the escape flux density (Hz -1 ) ratio of the Ly- 
man continuum (LyC, A ~ 900 A rest-frame) to the non- 
ionizing ultraviolet (UV, A ~ 1500 A rest-frame), 7?. eS c. 
This 7^. eS c is natura lly measured by direct observations (e.g., 
ISteidel et aljl200ll) . We compile all the observations of the 
LyC from galaxies and derive 7^ OS c at z < 3. On the other 
hand, 7Z CSC can be derived from the ionizing background 
intensity inferred from IGM absorption measurements via 
a cosmological radiative transfer model. We derive TZ CS c at 
z = 0-6 based on the recent reports of the background in- 
tensities. 

The cosmological parameters assumed in this Letter are 
H = 70.0 km s" 1 Mpc -1 , Q m = 0.3, and Q A = 0.7. 
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2 FORMULATION 



The definition of / csc in this Letter is 



/csc — 



LyC 



(1) 



LyC 



where J^yC an d ^Lyc are the intrinsic and escaping LyC flux 
densities, respectively. To obtain F{^ c , we have to correct 
the observed LyC flux density F£yQ for the photoelectric 
absorption by the neutral hydrogen remained in the IGM: 



riCSC 

^LyC 



F] 



LyC' 



; T LyC 



(2) 



We may estimate i^'ic fr° m the observed UV flux density 



pint 
^LyC 



lZ in tF vv e uv 



(3) 



where T^int = ^Lyc/^uv is the intrinsic LyC-to-UV flux den- 
sity ratio and Tyy is the UV dust opacity in the interstellar 
medium of a galaxy. Therefore, we have 



and 



F 



obs „ 

LyC r T IG " 
J £ LyC 



pobs 
UV 



(4) 



(5) 



This 7?. eS c is called the LyC-to-UV escape flux density ra- 
tio in this Letter. As found in equation (5), we can obtain 
7£esc from the observed LyC-to-UV flux density ratio with a 
correction for the IGM absorption ag ainst the LyC. This is 
similar to the approach introduced bv lSteidel et alJ i200ll) . 
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On the other hand, 1Z CSC can be obtained from the ob- 
served ionizing background intensity through a cosmological 
radiative transfer model. We describe the different aspects 
of this estimation in the four following sub-sections. 

2.1 Cosmological Lyman continuum transfer 

The mean specific intensity at the observed frequen cy v a as 
seen by an observer at redshift z is given by (e.g.. IPeeblesI 
1993) 



10.0 



J{v a , z ) = 



(1 + z f 
4ir 



dzf z {z)p{v,z)e-^°^\ (6) 



where the frequency v — v a (\ + z)/(l + z Q ), dl/dz is the 
line element, p is the emissivity per unit comoving volume, 
and T e s is the effective IGM opacity. Although the upper 
limit of the integral could be infinity, we set it to z\. The 
line element is 



dl , - c 
Tz^ Z) = H(z)(l + z) 



(7) 



where c is the speed of light and H is the Hubble parameter. 
The comoving emissivity can be written as 



p{v, z) = pqso 0, z) + p ga i {v,z), 



(8) 



where pqso and p ga i are the QSO and galactic emissivities, 
respectively. The IGM emissivity c an be omitted (see §2.3). 
The effective IGM opacity is (e.g-. lParesce et aljll98(i|) 



T e ff (^o, Zo, z) = / dz 



dNm 0N^0z^ 



(9) 



where d 2 N '/dNuidz is the cloud number on a line of sight 
per unit redshift z interval and per unit H I column density 
JVhi interval, and N\ and N u are the lower and upper limits 
of the column density of the IGM clouds. The cloud optical 
depth is r = <Jn(i'')Nm with cth the hydrogen cross-section 
at the frequency v 1 — v a (l + z')/(l + z D ). The frequency 
dependence of an for the LyC is assumed to be oc v~ A . 



2.2 IGM opacity 

The effective IGM opacity is determined by the cloud 

number function which can be expressed as (e.g., 
iMiralda-Escude fc Ostrikejll99(f) 



d 2 N 

dN m dz 



oc N H 



(10) 



If we assume a single index j3 (1 < /? < 2) for all A?hi 
and z, and Nicth <C 1 and N u aji ;$> 1, equation (9) can be 
approximated to ijZuo fc Phinnev|ll993l : llnoue et, al.l l2005) 



T e ff(^o, 



r(2-p)A[N lC m(v'w- \i+z'ydz' , (ii) 



where F is the usual Gamma function, A is the number 
of the IGM clouds with the column density between JVi 
and N u on a line of sight per unit redsh ift interval, and 
v' = v {l + z')/(l + z ). According to IWevmann et alJ 
lll998ft and lKim etaD ll200^ . we assume j3 = 1.5, (^,7) = 
(34.5,0.2) for 2 1.1 and (A, 7) = (6.3,2.5) for 2 > 1.1, 
and N\ = 4.4 x 10 13 cm -2 . Note that we have assumed the 
same number evolution along the redshift (the index 7) for 




Figure 1. IGM Lya opacity. The solid and dashed lines are the 
IGM Lya opacities calculated by equation (11) (see the text for 
the detail). The circles and squares with error-bars are the ob- 
served opacity taken from Songaila (2004) and Fan et al. (2006), 
respectively. 



the Lya forest and for the Lyman li mit syst e m (an d also for 
the damped Lya clouds), whereas iMadaul il995l) adopted 
di fferent indices for these clouds. As shown in the appendix 
of Ino ue et alJ 1120051) . the same 7 for all clouds is compatible 
with the recent obse rvations of the Lyman limit systems by 
iPeroux et ail J2003|) . 

While the IGM opacity model is used in this paper for 
the LyC, we display in Fig. 1 the opacity predictions at Lya 
line where comparisons with observations are po ssible. The 
Lya cross-section is based on Wi ese et alJ fiTi lil: : . Note that 
the same IGM clouds contribute to the LyC and Lya opac- 
ities. Our model (solid line) shows a very good agreement 
with the data at z < 5 but a smaller opacity than the data 
at z > 5. An explanatio n is that th e param eters of the cloud 
distribution taken from lKim et alJ <l2002l) who analysed the 
Lya forest at z < 4. If we assume (.4,7) = (1.3,3.5) for 
z > 4, the model (dashed line) shows a very good agree- 
ment with the data. 



2.3 Emissivities 



The LyC emissivity from galaxies can be expressed as 

Pgal(f, Z) = f(v)1Z CBC (z)pw(z) , (12) 



where puv is the observed galactic UV emissivity per unit 
comoving volume, and f(u) is the LyC frequency depen- 
dence of galaxies. We assume that f{v) oc v~ 2 (e.g., 
iFioc fc Rocca-Volmerangelll997l) . independently of the red- 
shift. The effect of the spectral slope on the estimated 1Z CSC 
is 15% at most if the index is changed from to —4. A 
redshift dependence of 7^. CS c is introduced. 

A key feature of equation (12) is that the LyC emis- 
sivity is directly related to the observed UV emissivity with 
the only parameter lZ esc that we wish to determine. In other 
words, our estimate for R eBC is free from uncertain dust cor- 
rection and from the intrinsic LyC production rate (or 7^mt) 
which depends on the initial mass function. 
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Figure 2. Observed total galactic UV emissivities per unit co- 
moving volume. The circles, asterisks, diamonds, triangles, and 
squares are taken from Schiminovich et al. (2005), Sawicki & 
Thompson (2006), Iwata et al. (2006), Ouchi et al. (2004), and 
Bouwens et al. (2006), respectively. The solid and dotted lines 
are the standard and the high-emissivity cases assumed in the 
calculation, respectively. 



The Pu*v values are obtained from the integration to 
zero luminosity of UV luminosity functions compiled from 
the literature. They are shown as a function of redshift in 
Fig. 2. The error-bars indicate observational 1-a errors. For 
practical purpose, we express p^y as the following functional 
form, Puv(^) = Pvvg(z), with 



9{z) 



. 2.7 , 



1+g 
. 4.0 



(0 < z < 1.7) 
(1.7 < z < 3.0) 
(3.0 < z) 



(13) 



The normalizati on adopted is p* TV = 2.7 x 10 2 6 erg s _1 Hz -1 
Mpc~ 3 based on lSawicki fc Thompson! J2006I) . Puv OI equa- 
tion (13) is displayed as the solid line in Fig. 2 and called 
the standard case. Because of the scatter at z > 3, we con- 
sider another case, the high-emissivity case, with g(z) — 1 
for 2^3, displayed in Fig. 2 as the dotted line. 

The QSO emissivity is given by the prescription in 
iBianchi et al] <l200l]) who obtained the LyC emissivity 
of QSOs b a sed o n th e optical l umino sity functions of 
IBovle et al.l feOOOh and IFan et al] ]200l]) and an average 
QSO spec trum. We here adopt a QSO spectrum parame- 
terized by Madau, Haar dt7fc Reesl lll999|) . 

The ionized IGM can contribute to a f raction of the LyC 
emiss ivity by the recombination process iHaardt fc Madaul 
Il996f). However, it is obviously secondary. Thus, we omit 
it. iMiniati et al] 11200411 suggested a significant contribution 
to the LyC by the thermal emission from hot gas shock- 
heated by the cosmological structure formation. In spite of 
its hard spectrum, however, the estimated hydrogen ioniza- 
tio n rate is an order o f magnitude smaller than that reported 
bv lBolton et al] <l2005|) at z = 2-4 an d is less than one-third 
of that reported by Fan et all feOOfil) at z = 5-6. Thus, we 
also omit the thermal emission by the structure formation. 
The omission of the IGM emissivity leads to an overestimate 
of IZ-csc that should be less than a factor of 2. 



2.4 Estimation of the LyC-to-UV escape flux 
density ratio 

If we enter the IGM opacity model described in §2.2 and the 
LyC emissivities described in §2.3 into equation (6), we can 
predict the background intensity at the Lyman limit as a 
function of lZ esc . This latter quantity can be evaluated by a 
comparison between the predicted and observed intensities. 

The observ ed intensities are o btained in the following 
ways. At z < 4, IScott et al.l i2002T) presented intensities at 
the Lyman limit based on their ob serva tions of th e QSO 
proximity effect. iBolton et ail J2005I) and lFan et"afl ]2006ft 
estimated hydrogen ionization rates (Fhi) at z = 2-4 and 
at z = 5-6, respectively, based on the observed IGM Lya 
opacity and a cosmological hydrodynamical simulation. The 
ionization rate Thi is 



hv 



dv 



47TCTL Jh 

h(a + 3) 



(14) 



where the subscript L means the quantity at the Ly- 
man limit, h is the Planck constant, and a is the 
power-law index when we assume a power-law back- 
ground radiation (J tx v~ a )- Assuming a — 2, we have 
Jl/W' 21 erg s" 1 cm" 2 Hz" 1 s r" 1 = r m /2 39 x 10" 12 s" 1 
with o l = 6.30 x 10" 18 cm 2 (IOsterbrocklll989l) . The un- 
certainty on Jl resulting from the adopted a is a fac- 
tor of (a + 3)/5. Wi th this conv e rsion we obtain Jl 
from the Thi values of lBolton et al] ]2005j) and IFan et al] 
feOOrjft . We should note here that IFan et al] J2006T) de- 
termined their absolute Thi values so as to fit those of 
iMcDonald fc Miralda-Esc udel (l20 0lh wh ich are somewhat 
smaller than those of lB'oltonet al] ]20o'5T) . 

Tiosc is determined by a comparison of these observed 
intensities at the Lyman limit (Jj^ s ) with the calculated 
theoretical intensities as 



tQSO 



QSO 



z)+Jt l (z,K e , 



(15) 



where and J£ are the Lyman limit intensities origi- 

nating from QSOs and galaxies, respectively. This compar- 
ison is performed by a non-parametric way from high to 
low redshift. Calculating J£ al (and then 7^. sc) at a red- 
shift, we use the TZ CSC values at redshifts larger than the 
redshift. For three different references of the observed inten- 
sities, we estimate 7?.osc independently. The assumed initial 

and 6.0 fo r the data 
l200Efl . and lFan et all 



reds hifts in equation ( 6 ) are 5.0, 5.0, 
fromlScott et al] j2002l) . lBoiton et al] 

(2006), respectively. This choice does not affect the results 
if we take an enough high redshift because the IGM opacity 
is very large at high redshift. In addition, we note that we 
assumed the higher opacity case at z > 4 (shown in F ig. 1 
as dashed line) only for the data from lFan et al] <|2006T) . 



3 RESULTS 

The values of 7Z eS c obtained at different redshifts either from 
direct observations or from indirect estimations (§2.4) are 
displayed in Fig. 3. 

Their comparison reveals an evolution of lZ CS c, with 
lZ ea c getting larger at higher redsh i fts. If we put confidence 
in the upper limit of lMalkan et al] i200St) at z ~ 1, which is 
an average value of 7?. e sc for 11 star-forming galaxies and is 
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Figure 3. Escape flux density ratio of the ionizing to the non- 
ionizing UV photons. The symbols and arrows with a label are 
measurements and upper limits from direct observations: 105 
(Inoue et al. 2005), S01 (Steidel et al. 2001), G02 (Giallongo 
et al. 2002), M03 (Malkan et al. 2003), L95/H97 (Leitherer et 
al. 1995, Hurwitz et al. 1997), D01 (Deharveng et al. 2001), and 
B06 (Bergvall et al. 2006). The circles and large downward arrows 
are the estimated value and upper limits based on the ionizing 
background intensities reported by Scott et al. (2002). Only upper 
limits are obtained at lower redshifts because the intensities from 
QSOs are almost equal to or exceed the observed intensities (see 
also Fig. 4). The triangles and squares are the estimated values 
based on the ionization rates reported by Bolton et al. (2005) and 
Fan et al. (2006), respectively; the standard emissivity of galaxies 
is assumed for the filled symbols, and the high-emissivity case is 
assumed for the open symbols. The dashed line is a possible fit 
to the evolution of 1Z CBC that will be used in Fig. 4. 



Figure 4. Lyman limit intensity of the background radiation. 
The circles and arrows at z ~ are measurements taken from 
Scott et al. (2002). The triangles and squares are converted from 
the ionization rates of Bolton et al. (2005) and Fan et al. (2006), 
respectively. The dashed line is the intensity from the QSOs 
(Bianchi et al. 2001). The solid and dotted lines are the total 
and galactic intensities, respectively. The standard galactic emis- 
sivity and an evolution of the escape flux density ratio shown in 
Fig. 3 as the dashed line are assumed. 



7l ea c shown in Fig. 3 as the dashed line, we obtain the dotted 
line for the galactic contribution and the solid line for the 
total intensity. The dashed line is the QSO contribution. 
The main contributor to the ionizing background radiation 
changes from galaxies to QSOs at z ~ 3. 



less affected than a single measurement by the randomness 
of the LyC escape, we can obtain a possible evolution, for 
example, as 

Tl csc oc (1 + z) 5 , (16) 

from z — 1.5 to z — 4.0, which is shown as the dashed 
line in Fig. 3. Since the uncertainties are still very large, 
this rather strong evolution is just an example. However, 
we note that an evolution of lZ asc is also found even in the 
high-emissivity c ase (open symbo ls) . Such an increasing 1Z CSC 
was suggested by Mciksin] (l2005t) in his Fig. 1 although the 
redshift coverage is smalle r than ours. 

The measurement bv lSteidel et al] (1200 it) is somewhat 
larger than our estimations. This may be caused by their 
sample selection bias. Indeed, their sample galaxies are taken 
from the bluest quartile in the UV colour among Lyman 
break galaxies. At lower redshif ts, we have only uppe r limits, 
except for the measurement bv lBergvall et af] <|2006T) . and a 
large dispersion. This is probably due to the small number of 
the observed g alaxies. Indeed, only 7 galaxies have been ob- 
served to date I Leither er et al ]|l995HDeharveng et al]l200l] 
Berg vall et aljl2006T) . The galaxy observed bv lBergvall et al] 
( 2006) (Haro 11) may be an exceptional one in the local Uni- 
verse as discussed in their paper. 

Fig. 4 shows a comparison of the calculated ionizing 
intensity with the observed ones. Assuming the evolution of 



4 SUMMARY AND DISCUSSION 

We have discussed the escape flux density ratio of the LyC 
to the non-ionizing UV (1Z CSC ) in a new framework for com- 
paring the direct observations of the LyC from galaxies with 
the cosmic ionizing background intensity. We have found a 
redshift evolution of TZ CEC , with larger 1Z CBC at higher red- 
shifts. Here, we translate !Z aS c into the escape fraction (/esc) 
of ionizing photons from galaxies. 

By equation (4), we can convert 7£ e sc into / csc if we 
have 72.i n t and Tjjv - For normal stellar populations with 
the standard Salpeter initial mass function and mas s range 
0.1-100 M , TZiut = 0.2-0.3 (e.g.. Ilnoue et aljl2005l) . From 
the observed UV slope, t'v" = 1-2 for the Lyman break 
galaxies at z ~ 3 (e.g., IShaplev et aljl2003l) although the 
uncertainty is very large. In this case, we have 7^-esc ~ /esc- 
If this relation is valid in the range z — 0-6, our result means 
that /csc increases from a value less than 0.01 at z < 1 to 
about 0.1 at z > 4. 

If we explain the evolving lZ aS c by a constant / CS c, we 
have two possibilities according to equation (4); (1) an evolv- 
ing TZ- ln t and (2) an evolving Txjy 1 . The evolving lZ ln t case 
needs an order of magnitude larger production rates of the 
LyC relative to the UV at z > 4 than at z < 1. This 
might be realized by a top-heavy stellar initial mass func- 
tion as suggested for zero metallicity, Population III (PopIII) 
stars. Very massive (> 100 Mq) PopIII stars show almost 
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black-bo dy spectra with an effective tempe rature about 10 5 
K (e.g., iBromm. Kudritzki. fc Loebll200lT) and would give 
7?.int — 2. However, such a transition of the mass function 
would occur at a much higher redshift than z ~ 2-4 be- 
cause it is thought to hap pen at a very low met allicity like 
~ 10~ 6 Solar value (e.g lschneider et al.ll200rj) . Although 
Ijimenez fc Haimanl j200d) proposed that 10-30% of stars in 
z = 3-4 star-forming galaxies are massive PopIII stars, we 
would need more such stars to obtain an enough large 7?4nt. 
The evolving r^/ 1 case needs an increment of 2-3 mag of 
from z < 1 to z > 4. Such a systematic increase of 
Jt^ 1 for galaxies at high redshifts is not reported, whereas 
iNoll et all (120041) reported an opposite trend for z = 2-4 (it 
may be caused by a selection bias in high redshift sample). 
We are therefore left with an evolution of / csc , with a larger 
/ esc at a higher redshift. 

The escape of the LyC is probably a random phe- 
nomenon. Thus, the evolving / csc suggests an increase of 
the fraction of galaxies showing a large escape. According 
to theoretical studies, a smaller galaxy can show a larger 
/ csc because of g alactic winds ijFuiita et alJ 1200*31) and/or 
champagne flows IjKitavama et al JI2004T I . Thus, the evolving 
/csc may suggest a decrease at higher redshift of the aver- 
age scale of galaxies, which is consistent with the cold dark 
matter scenario. The evolving / csc also suggests a morpho- 
logical evolution of galaxies; at higher redshifts, a large frac- 
tion of galaxies would show very disturbed structure caused 
by galactic winds and/or champagne flows, yielding a large 
escape of the LyC. 

We note two issues which should be assessed in future. 
If there are a large nu mber of t ype II QSOs at high red- 
shifts as suggested by iMeiksinl i2006l) . the required 1Z eS c 
may become much smaller. If the fraction of galaxies having 
a low-luminosity active galacti c nucleus increa ses towards 
high redshifts as suggested by IMeiksinl i2005l) , a scenario 
with an evolving 7?.mt without an evolving / es c may become 
possible. 

Finally, we strongly encourage new direct measurements 
of the LyC from galaxies to confirm/reject the evolving 7?.csc. 
For z ~ 3, the ground-based large telescopes are useful, and 
for z ~ 1, the GALEX satellite is appropriate. 
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